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Abstract 

We consider an application of modification of our variational-wavelet approach to some 
nonlinear collective model of beam/plasma physics: Vlasov/Boltzmann-like reduction from 
general BBGKY hierachy related to modeling of propagation of intense charged particle 
beams in high-intensity accelerators and transport systems. We use fast convergent multiscale 
variational-wavelet representations for solutions which allow to consider polynomial and ra- 
tional type of nonlinearities. The solutions are represented via the multiscale decomposition 
in nonlinear high-localized eigenmodes (waveletons). In contrast with different approaches 
we do not use perturbation technique or linearization procedures. 
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Abstract 

We consider an application of modification of our va- 
riational-wavelet approach to some nonlinear collective 
model of beam/plasma physics: Vlasov/Boltzmann-like re- 
duction from general BBGKY hierachy related to modeling 
of propagation of intense charged particle beams in high- 
intensity accelerators and transport systems. We use fast 
convergent multiscale variational-wavelet representations 
for solutions which allow to consider polynomial and ratio- 
nal type of nonlinearities. The solutions are represented via 
the multiscale decomposition in nonlinear high-localized 
eigenmodes (waveletons). In contrast with different ap- 
proaches we do not use perturbation technique or lineariza- 
tion procedures. 

1 INTRODUCTION 

We consider applications of numerical-analytical tech- 
nique based on modification of our variational-wavelet 
approach to nonlinear collective models of beam/plasma 
physics, e.g. some forms of Vlasov/Boltzmann-like reduc- 
tions from general BBGKY hierarchy (section 2). These 
equations are related to the modeling of propagation of 
intense charged particle beams in high-intensity accelera- 
tors and transport systems [1], [2]. In our approach we 
use fast convergent multiscale variational-wavelet repre- 
sentations, which allows to consider polynomial and ra- 
tional type of nonlinearities [3]-[16], [17]. The solutions 
are represented via the multiscale decomposition in non- 
linear high-localized eigenmodes (generalized Gluckstern 
modes, in some sense), which corresponds to the full mul- 
tiresolution expansion in all underlying hidden time/space 
or phase space scales. In contrast with different approaches 
we don't use perturbation technique or linearization pro- 
cedures. In section 3 after formulation of key points we 
consider another variational approach based on ideas of 
para-products and nonlinear approximation in multiresolu- 
tion approach, which gives the possibility for computations 
in each scale separately [18]. We consider representation 
(4) below, where each term corresponds to the contribution 
from the scale i in the full underlying multiresolution de- 
composition as multiscale generalization of old (nonlinear) 
SF approach [1]. As a result, fast scalar/parallel modeling 
demonstrates appearance of high-localized coherent struc- 
tures (waveletons) and pattern formation in systems with 
complex collective behaviour. 



2 VLASOV/BOLTZMANN-LIKE 
REDUCTIONS 

Let M be the phase space of ensemble of N particles 
(dmxM = 6N) with coordinates xi = (qi,Pi), i = 
1,...,N, Qi = (qlqlcft) £ R 3 , Pi = (pj,p?,p?) G 
R 3 with distribution function Dn(xi, . . . , xn] t) and 

F N (xi, . . .,x N ;t) = }^D N (xi, . . . ,x N ;t) (1) 



be the N-particle distribution functions (Sn is permutation 
group of N elements). For s=l,2 we have from general 
BBGKY hierarchy: 

dFxjx^t) pi 9 

oT 1 ■^—Fi(xi;t) (2) 

ot to oqi 

= - ( dx 2 L 12 F 2 (x!,x 2 ;t) 
v J 

dF 2 (xi,x 2 ;t) (pi d p 2 d \ 

r I — "a 1 a Ll2 I V) 

ot V to oqi m oq 2 / 

F 2 (xx,x 2 ;t) = - [ dx 3 (L 13 + L^F^xx^x^t) 
v J 

where partial Liouvillean operators are described in [17]. 
We are interested in the cases when 

k 

F k (xi, . . .,x k ;t) = Y[Fi(xi;t) + Gk(x%, . . - ,x k ;t), 



where G k are correlation patterns, really have additional 
reductions as in case of Vlasov-like systems. Then we have 
in (2), (3) polynomial type of nonlinearities (more exactly, 
multilinearities). 

3 MULTISCALE ANALYSIS 

Our goal is the demonstration of advantages of the fol- 
lowing representation 



F = J2 61 F 



(4) 



iez 
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for the full exact solution for the systems related to equa- 
tions (2), (3). It is possible to consider (4) as multiscale 
generalization of old (nonlinear) SF approach [1]. In (4) 
each 8 l F term corresponds to the contribution from the 
scale i in the full underlying multiresolution decomposi- 
tion 

. . . C V-x C V C Vi C V 2 . . . (5) 
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of the proper function space (L 2 , Hilbert, Sobolev, etc) to 
which F is really belong. It should be noted that (4) doesn't 
based neither on perturbations nor on linearization proce- 
dures. Although usually physicists, who prefered computer 
modelling as a main tool of understanding of physical re- 
ality, don't think about underlying functional spaces, but 
many concrete features of complicated complex dynamics 
are really related not only to concrete form/class of opera- 
tors/equations but also depend on the proper choice of func- 
tion spaces where operators actully act. Moreover, we have 
for arbitrary N in the finite N-mode approximation 

JV 

pN = Yl 6%F (6) 

i=l 

the following more useful decompositions: 



{F(t)}= W S or {F(t)} = V @W j , (7) 

— oo<j<oo j=0 

in case when Vq is the coarsest scale of resolution and 
where Vj+i = Vj ® Wj and bases in scale spaces Wi(Vj) 
are generated from base functions tp(ip) by action of affine 
group of translations and dilations (the so called "wavelet 
microscope"). The following constructions based on varia- 
tional approach provide the best possible fast convergence 
properties in the sense of combined norm 

\\F N+1 -F N \\<e (8) 

introduced in [17]. Our five basic points after functional 
space choice are: 

1 . Ansatz-oriented choice of the (multidimensional) ba- 
ses related to some polynomial tensor algebra. Some 
example related to general BBGKY hierarchy is con- 
sidered in [17]. 

2. The choice of proper variational principle. A few 
projection/ Galerkin-like principles for constructing 
(weak) solutions are considered in [3] - [16]. It 
should be noted advantages of formulations related to 
biorthogonal (wavelet) decomposition. 

3. The choice of bases functions in scale spaces Wj 
from wavelet zoo. They correspond to high-localized 
(nonlinear) oscillations/excitations, coherent (nonlin- 
ear) resonances, etc. Besides fast convergence prop- 
erties of the corresponding variational-wavelet expan- 
sions it should be noted minimal complexity of all un- 
derlying calculations, especially in case of choice of 
wavelet packets which minimize Shannon entropy. 

4. Operators representations providing maximum sparse 
representations for arbitrary (pseudo) differential/ in- 
tegral operators df/dx, d n f/dx n ,J T(x, y)f(y)dy), 
etc [17]. 

5. (Multi)linearization. Besides variation approach we 
consider now a different method to deal with (polyno- 
mial) nonlinearities. 



We modify the scheme of our variational approach in such 
a way in which we consider different scales of multiresolu- 
tion decomposition (5) separately. For this reason we need 
to compute errors of approximations. The main problems 
come of course from nonlinear (polynomial) terms. We 
follow according to the multilinearization (in case below - 
bilinearization) approach of Beylkin, Meyer etc from [18]. 
Let Pj be projection operators on the subspaces Vj (5): 

(^■/)W-E < /'^ > ^,feW (9) 

k 

and Qj are projection operators on the subspaces Wj\ 
Qj = Pj-i — Pj. So, for u £ L 2 (R) we have Uj = PjU 
and Uj £ Vj. It is obviously that we can represent u 2 in the 
following form: 

n n 

u 2 = 2j2(Pju)(Qju) + J2(Qju)(Qju) + u 2 n (10) 

2=1 2=1 

In this formula there is no interaction between different 
scales. We may consider each term of ([!(]) as a bilinear 
mappings: 

Kw ■■ v j x W J -» ^ 2 ( R ) = VjtBj^jWj, (11) 

■ Wj x Wj -» L 2 (TL) = Vj ® f >j Wy (12) 

For numerical purposes we need formula ( |lO| ) with a finite 
number of scales, but when we consider limits j ^ oo we 
have 

u 2 = ^2(2P jU + Qj U )(Qju), (13) 
2'ez 

which is para-product of Bony, Coifman and Meyer [18]. 
Now we need to expand ( |lC)| ) into the wavelet bases. To 
expand each term in ( |Io| ) we need to consider the integrals 
of the products of the basis functions (7), e.g. 

/oo 
^{{xW k ,{x)4{x)dx, (14) 
-oo 

where j' > j and 

ip{(x) = 2-^ 2 i){2~ j x - k) (15) 

are the basis functions (7). For compactly supported 
wavelets 

M$ ww (k,k'J) = for |fe-fe'|>A;o, (16) 

where ko depends on the overlap of the supports of the basis 
functions and 

W www (k - k', 2 r k -1)\<C- 2- rXM (17) 

Let us define jo as tne distance between scales such that for 
a given e all the coefficients in ( |l7| ) with labels r = j — j', 
r > jo have absolute values less than e. For the purposes 
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Figure 1: N = 1 waveleton contribution to (6). 




Figure 2: Stable pattern. 



of computing with accuracy e we replace the mappings in 



M^ w : W s x W 3 



V, 



3 *&3<j'<j( 



V, 



3 ^j<j'<30 W f 



(18) 
(19) 



Since V 3 ®j< r < jo w f = V jo ^, V 3 c v^, m C 
Vj -x we may consider bilinear mappings ( |l8[) , (Uy) on 
Vj -i x Vj -i. For the evaluation of (p"8[), ( ]l9[ ) as map- 
pings Vj -i x V5 -i — > Vj _i we need significantly fewer 
coefficients than for mappings (|l8]), (|T^). It is enough to 
consider only coefficients 



M(k, k', I) = 2- j /' 2 / tp{x - k)<p(x - k')cp(a 



where ip(x) is scale function. Also we have 

M(k,k',i) = 2- j / 2 M (k-£,k' -I), 

where 



£)dx, 
(20) 

(21) 



M (p, g) 



(^(x — p)f{x — q)ip(x)dx (22) 



Alo(p,q) satisfy the standard system of linear equations 
and after its solution we can recover all bilinear quanti- 
ties (14). Then we apply some variation approach from 
[3]-[16], but in each scale separately. So, after applica- 
tion of points 1-5 above, we arrive to explicit numerical- 
analytical realization of representations (4) or (6). Fig.l 



demonstrates the first contribution to the full solution (6) 
while Fig. 2 presents (stable) pattern as solution of system 
(2)-(3). We evaluate accuracy of calculations according to 
norm introduced in [17]. 
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